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Abstract. We investigate the spectrum of a quasi-monoen- 
ergetic ensemble of relativistic electrons, especially for the 
mildly relativistic case, and discuss the effect of inclination of 
the magnetic field on the emissivity. We apply the exact theo- 
retical description to the spectrum of the radio source Sgr A* 
which is located at or very close to the dynamical center of 
our Galaxy. We find that the radio-MIR spectrum can be repro- 
duced well, but that the resulting self-comptonized X-ray flux 
is much smaller than the observed one. 



Key words: Galaxy: center - Radiation mechanisms: non-ther- 
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1. Introduction 

Observations of the galactic center radio source Sgr A* in the 
last few years (e.g. Zylka et al. 1995) have raised the inter- 
est in detailed investigations of synchrotron spectra produced 
by mildly relativistic elec trons (e.g. Melia 1994, Mahadevan et 
al. 1996). [ Melia (1994) ] proposed synchrotron radiation from 
thermal electrons in a radial inflow but failed to explain the lat- 
est sub-mm data. Narayan et al. (1995) presented a model for 
the entire spectral range from radio to X-ray frequencies based 
on advection dominated accretion. 

In our view, only the radio observations can be taken as 
firm detections and the best fit from accretion models with LTE 
is based on a isothermal sphere emitting optically thin syn- 
chrotron radiation. In a previous paper (Beckert et al. 1996, 
hereafter referred to as B96) we showed that a homogeneous 
blob of relativistic electrons penetrated by a magnetic field is 
sufficient to explain the observed spectrum. 

In the present paper we describe the calculation of syn- 
chrotron spectra of electrons with various distributions in mo- 
mentum space as seen by an arbitrary observer. The basic equa- 
tions are presented in Sect. p| and A ppendix A. We apply the 
general formulae (e.g. [Bekefi 1966]) without approximations 
to a source of mildly to highly relativistic electrons. In Sect. |] 



we discuss in detail the synchrotron spectrum of monoenergetic 
electrons and compare the results with the standard treatment 
for highly relativistic electrons. In Sect. Qwe provide a criterion 
for quasi-monoenergetic distributions and show the deviations 
from a simple approximation to the synchrotron spectra which 
is used to identify optically thin monoenergetic synchrotron 
sources (Duschl & Lesch 1994, Reuter & Lesch 1996). In the 
subsequent Sect. [5] we derive the synchrotron flux of homoge- 
neous spherical sources with optical thick-thin transition. Sect. 
^ deals with the application of quasi-monoenergetic electron 
distributions to the observations of Sgr A*. We interpret the 
low-frequency turn-over (B96) in terms of synchrotron self- 
absorption (SSA) and determine the physical source parame- 
ters under the assumption of energy equipartition between elec- 
trons and magnetic field. An unavoidable consequence of the 
synchrotron model is the self-comptonization of the emitted 
radiation (SSC) which emerges as X-rays, briefly discussed in 
Sect. 0. We summarize our results in Sect. §. 

2. Synchrotron Theory 

The theory of synchrotron emission has been reviewed several 



times ([Bekefi 1966 1, Ginzburg & Syrovatsky 1969). Here we 



briefly discuss the basic formulae for the coefficient of spon- 
taneous emission and for the observed power of synchrotron 
radiation of single electrons. The fundamental work was done 
by Schwinger (1949) and Westfold (1959), but later on the dif- 
ference in observed and emitted power of single electrons was 
recognized by Epstein & Feldman (1967) and Scheuer (1968). 
We use a formulation of synchrotron emision based on the re- 
tarded potentials to calculate the recieved power seen by a dis- 
tant observer. Radiation losses are neglected and we assume 
that binary encounters do not take place. The relativistic elec- 
trons are characterised by their energy E 
frequency 

qB 1 



7m 



,c , the orbital 
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in a locally homogeneous magnetic field B and a pitch angle 
a between the magnetic field and the momentum of the spi- 
raling electron. The electric and magnetic radiation field which 
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determine the flux density are peroidic in the retarded time of 
emission r 



t 



(2) 



Rather than performing a Fourier transformation to get the 
spectral emissivity of a single particle in a homogeneous mag- 
netic field we expand the radiation field in a Fourier series. The 
recieved power per steradian in a direction 9 to the magnetic 
field is derived in Appendix A and reads 

J2 00 



P = 



— 

2ttc ^ 



with a complex field vector 

muj s ( ^ sin a (mi/;) 

g m = — «(/3cosa — cos#)sin6? _1 J m (mij)) 

« V o 

depending on the geometric quantities 

sin a sin 9 



£ = 1 — j3 cos a cos 9 



(3) 
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The coefficient of spontaneous emission r\ v which is related to 
the emissivity by j v = J d 3 p £ f{p)rj v (p) is obtained with 
a resonance condition 5 (y — mVs(T), since we are dealing 
with emission at discrete frequencies v = mvg^ 1 . 
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'/i = I —, : I - (7) 

p sin 9 

This form of the spontaneous emission differs from the stan- 
dard form (e.g. [Bekefi 1966]) by an additional factor In 
the standard calculation a square of <5-functions appears. From 
our point of view, one of these is replaced by ojJ 1 and di- 
vided by T — to get the mean power from the total 
energy. But the frequncy v s — T^ 1 is not the observed fre- 
quency v — £i/ s of the radiation. This additional £ is contained 
in our Eq. (||). We will in see Sec. [| that this difference is offset 
by the normalization of the distribution function in calculating 
the spectral emissivity. 



3. Synchrotron spectra of monoenergetic electrons 

We describe the radiating electrons by a stationary distribu- 
tion function f e (a,p) in momentum space in a homogeneous 
magnetic field. The distribution in phase space is assumed to 
be separable and the complete distribution function looks like 
n(r)f e (a, p) with the volume density n(r). In order to recover 
the standard description of emissivity we must consider the dis- 
tribution / e in the observers frame. We take f e to be normal- 
ized so that the density of electrons in the source is contained 
in n(r). In principle, the distribution is a density and can be 
time dependent. While the distribution describes the electrons 
in the source at time r, we have changed the time variable to 
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Fig. 1. The spectrum of mildly relativistic electrons with 7 = 5 
seen by observers at 9 = 5°, 20°, 60°, 82°. We plot the re- 
ceived power vj v for a normalized distribution. The spectra ob- 
tain their maxima at higher frequencies as observers approach 
9 = 90°. The thickest solid line is the average of all observers 
corresponding to an isotropic distribution of magnetic field di- 
rections. 
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Fig. 2. The flux density distribution of Sgr A* fitted with mo- 
noenergetic electrons of 140 MeV corresponding to 7 = 275. 
The observed spectrum is taken from B96 where symbols with 
circles at both ends indicate variability. The solid curve cor- 
responds to the mean observer position and broken lines refer 
to 6 = 89°, 45°, 20°, 2° from high to low flux densities at 10 
GHz. For physical parameters, see model I in Tab.[[J 



the time of observation t. This transformation is reflected in 
the differential dt = £ dr and the distribution seen by the ob- 
server is therefore £f e (a,p). A more extended discussion of 
this argument is given in Appendix B. The integrated spectral 
emissivity in the observers frame follows from the power of a 
single electron stated above 
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The spectral emissivity is obtained from the spontaneous emis- 
sion given in Eq. (||) by evaluating the a-integration. The 5- 
function fixes a for each combination of [v, 9 ^ -| , to] unam- 
biguously by 

cosa m = — ^- . (9) 

pcost' 

Since a rn 6 [0, it] has to be a real number, we get a condition 
for the index m 

— (1 -P\cos0\) < to < — (1 + f3\cos9\) . (10) 
v s v s 

It turns out that the problem of deducing the spectrum of a dis- 
tribution of electrons in momentum space is reduced to one 
integral and a limited sum of Fourier coefficients : 



j u (Ct) = n(r) 
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The abbreviation g\ appeared previously in Eq. (Q). This form 
is an exact result with a m determined by Eq. (^). The sum can 
be evaluated directly for small frequencies v ~ v s and small 
momenta where only a few cyclotron lines emerge in the spec- 
trum. For frequencies where the range of indicies m for the 
harmonics contributing to the sum is wide, it is reasonable to 
transform the sum to an integral and evaluate it numerically. 

For a monoenergetic and isotropic distribution, the distri- 
bution function is fe{oi m ,p) = (47rpo) _1 ~ Po) an d the 
emissivity in Eq. ( p"l| ) consists of contributions from electrons 
described by (a m ,po) whose TO-th harmonic is shifted to the 
frequency v in the direction of the observer. The resulting spec- 
trum is shown for four different observers in Fig.|lJ for electrons 
with a common Lorentz factor of 7 = 5 corresponding to 2.56 
MeV. We see that the average of all observers is similiar to the 
emission in the direction 9 = 60° shown by the dashed line 
in the Figure. The received flux changes drastically when the 
direction to the observer becomes parallel to the magnetic field 
lines. The emission for 9 — 5° consists of overlaping low m- 
harmonics of electrons with small pitch angels. In addition, it 
demonstrates the transition from overlaping cyclotron lines to 
a continuous synchrotron spectrum. The flux measured by the 
observer is increased relative to the emitted power of the spi- 
raling electron which mainly radiates in the direction parallel 
to its momentum. Ginzburg et al. (1968) recognized that there 
also is a change in the total energy of the radiation field inside a 
fixed sphere in the observers frame containing the electron and 
the observer at its surface. 

The synchrotron spectra of monoenergetic electrons pro- 
vide a reasonable fit to the Sgr A* spectrum except for the turn- 
over below 1 Ghz. We have taken the average of all observer 
positions and a magnetic field strength of 10 G, comparable 
to the first model using monoenergetic electrons by Duschl & 
Lesch (1994). The ^-dependency is increased for high ener- 
gies as seen in Fig.|| The parameters are combined in Tab.|l]as 
model I and discussed in greater detail in Sect. ||. 

So far we have a description for the spontaneous emission. 
To get a complete spectrum, we must include induced emission 



and absorption inside the source. The cross section for syn- 
chrotron self-absorption (SSA) is basically given by 

r 2 r°° 

^=g^3 J dpp 2 (f( Pl )-f(p))e u (p) (12) 

where p\ is the momentum of the electron before absorption. 
The relativistic limit - which we assume for SSA - is simply 
P\c — pc — hv. The emissivity e u (p) is the ^-average of j v per 
electron. In principle the emissivity has to be split up into the 
two directions of polarisation present in Eq. (Q) and the cross- 
section must be calculated separately for both directions. Since 
we assume that the distribution of magnetic field directions is 
isotropic on a length scale much shorter than the source size, 
this distinction can be omitted without indroducing relevant er- 
rors. Using a Taylor-expansion for f(pi) — f(p) we get an ex- 
pression which closely resembles the standard formula and can 
easily be evaluated numerically. 

°v = o C , 3 / dp(px-p)^-f V — ^ , \g m \ 2 -(13) 

Sithv 6 J op ' mv s |cos0| 

All low frequency turn-overs seen in the spectra of Fig. |o]-[To| 
are due to sychrotron self-absorption. 



4. A criterion for quasi-monoenergetic distributions 

Untill now we have not defined what we call a quasi-mono- 
energetic electron distribution. To do this we compare the ex- 
act formulae given in Sec. || with a simple approximation for 
the emissivity of monoenergetic electrons in the energy range 
7 = 100 . . . 1000. This will allow us to derive analytic expres- 
sions for the emissivity of truncated power-law distributions. 
From these we get a criterion what a quasi-monoenergetic dis- 
tribution will be when synchrotron emission is considered. We 
approximate the synchrotron spectrum of a single electron by 

P0, 7 ) = P u 1 / 3 eiq>[-v/v c } (14) 

which is used in the interpretation of the radio spectra of Sgr 
A* (Duschl & Lesch 1 994) and the core of M8 1 (NGC303 1 ) by 
Reuter & Lesch (1996). The total luminosity of a single elec- 
tron is 

L = L B 2 {-fP sin a) 2 

with pitch angle a in the frame of the source and 
L = 1.5870 lO-^ergs- 1 

if the magnetic field is measured in Gauss. In the observers 
frame the shift of time intervals dt = £dr introduces a factor 
£ -1 as discussed in Sec. || and in Appendix B. For sufficently 
high electron energies and frequencies much larger than the 
orbital frequency, the radiation is beamed in the direction a w 
9 and we approximate £ w sin 2 a. The observed luminosity 
becomes independent of the pitch angle : 

L = L a B 2 j 2 

The scaling of the spectral power P{v, 7) is obtained from the 
luminosity and gives 

nBq^ 2 



i/ 7 



27TTOC 
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As a factor in the critical frequency v c we introduced a coef- 
ficient k which is usually taken as k = 1.5. We would like to 
compare the approximation with the exactly evaluated spectra 
and minimize the relative error. This suggests a smaller value 
for k ps 1.2 so that the error is less than 20% in the range 
100 < 7 < 1000 as shown in Fig. |3|. We restrict the error 
box in Fig. [| to frequencies around the maximum of the flux. 
At high frequencies v > hv c the spectrum drops nearly expo- 
nentially and the flux has decreased by two orders of magni- 
tude. At low frequencies the spectrum is changed by the pres- 
ence of cyclotron lines which provide a significantly smaller 
flux. But the absolute flux is small at low frequencies due to 
the approximated P u ~ v 1 ! 3 behavior and is often suppressed 
by self-absorption. In order to get analytic expressions which 
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Fig. 3. Relative error of the simple approximation Eq.(|T4|) com- 
pared with the monoenergetic case averaged over the observer 
positions 8. The crosses indicate the error for 7 = 1000 and the 
solid line for 7 = 100. The deviations at small frequencies are 
dominated by the discrete cyclotron lines present in the strict 
monoenergetic case. The magnetic field is taken as B = 10 G. 



determine what we want to call a quasi-monoenergetic elec- 
tron distribution, we investigate restricted power-laws /o7 _a 
for 71 < 7 < 72 and set the distribution to zero for other en- 
ergies. The width \ of the distribution is defined as \ — 72/71 
and is the main parameter in our investigation. The mean lumi- 
nosity is then 

Q 2 (2) 



(L) = n(r)L a B ( 7 )' 



Q{x) = 



0(1) <9(3) 



We normalize the distribution function so that n(r) gives the 
spatial density of the relativistic electrons. 

n(r,7) = n(r)/o7 _Q 



/o=7r 1 Q(l) (7) =7i 



m 

Q(2) 



We assume equipartition between magnetic and kinetic energy 
density. Despite the presence of a thermal plasma, we only con- 
sider the relativistic electrons in the kinetic energy. The spatial 
density 

n(r) = 2 ( 7 ) _1 

airmc 

is thus unambiguously determined by the magnetic field and 
the mean electron energy. In the discussion of quasi-monoen- 
ergetic distributions, we want to redefine the critical frequency 
given in Eq.(|l5|) by replacing the square of the relativistic 7- 
factor by the square of the mean (7 ) 2 

_ KBq(j) 2 



2-n:mc 
The mean emissivity is 

J,=J S 8/ V/ 3 ( 7 )-*G 
with 

3L 
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d7 7 -Q ~ 2/3 exp[-t7t/ /7 2 
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(17) 



If v < v c we find G w Q 5 / 3 (l) Q- 2 / 3 (2) Q-^l/S). To 
see the effect in the frequency range v rj (7)^0 we must 
evaluate the correct integral U.We introduce a new variable 



U = 



7 

v 
3Bq 



2 an find : 



ds s a / 2 " 7 /6 eX p[-s] 



Attitic 

The integral U can be expressed as the difference of two in- 
complete Gamma-functions 7(0, x) 



U 



V v 



a 1 

71 7t-e, 



1 



' 1 2 6^071 



2 6V o7l 2 

or as the difference between two confluent hypergeometric 
functions M(a, b, z) (Kummer functions) resulting in a com- 
pact form for G 

' - 1/3 - Q Q(72) 



G = 2e-Z Q 5 / 3 (l)Q- 2 / 3 (2) 
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Q( 7 ) = exp -v/vq (7 
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The changes in the synchrotron spectra become significant 
when the width of the distribution grows beyond a limit of say 
X = 5. In that case the effect of the broader distribution ex- 
ceeds the error introduced in the approximation Eq. (|l4j). The 
effect of the width on the approximative spectra are shown in 
Fig. H were the solid line correspond to x — 7. 
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For a homogeneous, spherical, optically thin source of ra- 
dius R we find for the received flux at a distance d 
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(18) 



Si/ 



with the coefficient n taken to be 1.2 . For completness we 
present the absorption coefficient for synchrotron self-absorp- 
tion in the same approximation as above for the emissivity. In- 
serting the mean emissivity from Eq. (|l6|) in Eq. ( |l2| ) we get 

a v = aoBiv-i ( 7 )"* Q 8 / 3 (l) £T 5 / 3 (2) J 
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We will not discuss synchrotron self-absorption further in this 
approximation. 



Fig. 4. Comparision of the approximation Eq. ( |l7[ ) to the nu- 
merical evaluation of Eq. (^lj) for a truncated power-law of 
width x — 3, a = 2 averaged over all observers. The approxi- 
mation is shown as a dotted line and the numerical equivalents 
are plotted for B = 10 G and 7 = 1000 (solid) and 7 = 100 
(dashed). The lower panel presents the relative error of the ap- 
proximation. The dotted lines mark the 0% and 10% level . 
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Fig. 5. Deviations of the synchrotron spectra of quasi- 
monoenergetic electrons from the strictly monoenergetic case. 
We took power-laws with a — 2 and plot G — exp[— vjv^ 
as defined in Eq. (p^). The deviations grow with \ = 
1.5,3,5,7,10,20. 



5. Aspects of radiative transfer in a homogeneous and 
spherical source 

We assume that the strengh of the magnetic field B, the den- 
sity of relativistic electrons n(r) and the energy distribution of 
the electrons are constant throughout the source. As a conse- 
quence, the emissivity e„ and the absorption coefficient a v do 
not depend on the location inside the source. This simplifies the 
radiative transfer so that we can get analytic expressions for the 
synchrotron flux at the surface. At first we find for the intensity 



T = io_ a Q(2) (7) U n 
" a (2 + a) Q(l) 7l J 



exp | 



*„s]) 



if s is the length of the path through the blob. Here we have 
again assumed energy equipartition as in Sec. If we integrate 
over the surface in the sky we find a total flux 
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(19) 



In the optically thick regime the flux is proportional to the area 
2irR 2 and in the optically thin case we perform a Taylor expan- 
sion and recover the simple form 



F„ 



Air 



-R 6 J V 1 
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to first order in in the optical depth Ra v . 
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Table 1. Physical parameters for the synchrotron source Sgr A* 
Model (E[MeV]) B [G] R [10 13 cm] n [10 3 cm -3 ] 



I 


140 


10 


1.25 


18.0 


11 


217 


2.7 


4.92 


0.84 


111 


222 


3.0 


4.86 


1.01 


IV 


210 


4.0 


3.75 


1.90 


V 


202 


4.0 


3.50 


1.97 



6. Application to the galactic center source Sgr A* 

The idea of monoenergetic relativistic electrons producing op- 
tically thin synchrotron radiation has been applied to Sgr A* by 
Duschl & Lesch (1994). Later on, Narayan, Yi & Mahadevan 
(1995) showed that thermal electrons in an advection domi- 
nated accretion flow can lead to quite similiar spectra, but an 
isothermal sphere seems to fit the data better. In B96 we dis- 
cussed power-laws for isotropic electron distributions of lim- 
ited width with infinitely steep cut-offs and emphasized the 
importance of the absorption mechanism for determining the 
source parameters. Here, we extend the distributions from the 
strict monoenergetic case discussed above to include more real- 
istic power-laws with different cut-offs and relativistic thermal 
distributions. The spectra in B96 were derived for a position 
angle of 8 = 60° and in Sect. [| we showed that this is a good 
approximation to the average of all possible position angles. 
Here we always take the average of all directions of magentic 
field lines. For a complete description of the source, we as- 
sume energy equipartition in order to relate the magnetic field 
strength with the number density of electrons. We calculate the 
spectral flux density according to Eq. ( p^ ) and adopt a distance 
of 8.5 kpc to Sgr A*. The model V in Tab. |l] is the best fit with 
a truncated power-law distribution as discussed in Sec. [| and 
B96. 



6.1. Power-laws with exponential cut-offs 

In B96 we introduced infinitely steep cut-offs for power-law 
distributions to achieve an easier understanding of the char- 
acteristics of synchrotron emission from quasi-monoenergetic 
electron distributions. They are not expected from acceleration 
mechanisms for electrons in a relativistic plasma. However, 
distribution functions with an exponential cut-off at high en- 
ergies and a rising power-law at low energies are expected for 
equilibrium states (e.g. Schlickeiser 1984). If we think of an in- 
jection mechanism at high energies like magnetic reconnection, 
we assume exponential decays at both ends of the distribution 
and a power-law in between. The distribution, normalized for a 
single electron, becomes 

_ _ Pi 

P . 
IPx 
Pi, 

The normalization is expressed in terms of modified Bessel 
functions K z (x) of order z. In contrast with power-laws cov- 
ering a wide range of electron energies, the spectral power a 



f(p) =K-Vexp 



K = 2( PlP2 ) — K 1+a 




(20) 



does not determine the emission spectrum and is assumed to be 
cr = —1 in our case. The width \ = P2/P1 is defined in the 
same manner as in Sec. f|. The best possible model obtained 
with a single homogeneous blob of such electrons is shown in 
Fig. H The physical parameters are given as model II in Tab. 
[l]. In contrast to the spectra with infinitely steep cut-offs, this 
model cannot account for the measured flux densities in the 
range 200 — 1000GHz. The power a and the width (here x = 3) 
are unimportant for this part of the spectrum. It is dominated by 
an exponential tail beyond p 2 . 



6.2. Thermal distributions 




Fig. 6. Synchrotron spectrum of electron distributions with ex- 
ponential high energy tails. Both thermal electrons of T = 
6 10 n K (solid line ; model III in Tab.|l]) and power-laws with 
exponential cutoffs (dashed line; model II) at high and low en- 
ergies can hardly be distinguished. 



Another closely related spectrum is produced by relativistic 
thermal electrons described by a distribution function 



f(p) = Kp 2 exp 



7 



as seen in Mahadevan et al. (1996) with the corresponding nor- 
malisation K to the number of electrons. The synchrotron spec- 
trum beyond the maxima at l.Z^Vg is dominated by the expo- 
nential decay in the high energy band of the electron spectrum. 
This is also the case for the distributions with exponential cut- 
offs discussed above. Obviously, both spatially homogeneous 
models suffer from the same disease, since they can not ex- 
plain the sub-mm data seen in Fig||. The parameters for the 
thermal distribution are listed as model III in Tab. [jj. 

6.3. A Gauss-Profile for the electron momenta 

In the previous example we saw that an exponential high en- 
ergy tail of the distribution is in contradiction to the assumption 
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of one single homogeneous blob of optically thin synchrotron 
emission. A steeper tail is present in a Gaussian profile 

f( p )=K-VeKp\-^^] (21) 
with 

for the norm of the distribution. We have included a phase space 



Table 2. Two component model for the Synchrotron source Sgr 

A* 



1.00 



0.10 



0.01 

10 




10 10 10 

v [GHz] 



Fig. 7. Synchrotron flux from electrons with a gaussian profil in 
momenta according to model IV in Tab. |] The width is varied 
in the range \ = 0.2 (dashed), 1.0 (solid), 1.92 (dotted) and 
produce spectra with increasing flux densities at 2 10 4 GHz. 

factor of p 2 which does not affect the resulting spectra very 
much. Again, the high energy tail is the most important part 
necessary to explain the sub-mm measurements and yet stay 
below the upper limits in the IR. We define the width of the 
Gaussian distribution \ as the FWHM Ap of the pure gauss- 
profile relative to the mean momentum p\ so that \ = ^p/pi 
. The momentum p2 is now determined by the width pi — 
0.6 xPi- Values greater than 2 are not reasonable for this def- 
inition of the width. The obtained fit corresponds to model IV 
of Tab.|. We see in Fig. ^ that this distribution is sufficient to 
explain the radio spectrum of Sgr A* for a wide range of widths 
uptox = 1-92. 



6.4. A core-shell model for Sgr A* 

Since distributions with an exponential high energy tail fail to 
explain the full spectrum, using the assumption of isotropy and 
homogeneity of the source, we can build up a two component 
source. There is good observational evidence for the existence 
of a central black hole of mass M w 2.5 10 6 A/ Q in the galactic 
center (Eckart & Genzel 1996) which powers Sgr A*. Accre- 
tion into this object implies a special radial structure suggested 
by Narayan, Yi & Mahadevan (1995). Our first step towards a 
radial structure is a core-shell model. 



Comp. (g[MeV]) B [G] 



R [10 13 cm] 



n [10 4 cm -3 ] 



1 

II 



155 
41.4 



2.0 
70 



5.5 
0.13 



0.064 
300 



While the observational situation is far from clear yet, al- 
most simultaneous observations of the spectrum of Sgr A* 
from cm to sub mm wavelengths (Falcke et al. 1997) seem to 
indicate the possibility of some execess submm flux in com- 
parison to the flux predicted by the homogeneous model. 

Our core-shell model consists of an optically thin extended 
source and a highly self-absorbed compact component. This 
is not the result of a hydrodynamical calculation, but the best 
description to explain the observations. The extended part, 



1.00 r 



0.10 r 



0.01 L 

10* 




Fig. 8. The spectrum of a core-shell model of Sgr A* consisting 
of two thermal emitting regions which are homogeneous in it- 
self. The model parameters are given in Tab.|| for the extended 
shell as part I and the self-absorbed compact core II. Both com- 
ponents are drawn separately as dotted lines. 

component I, dominates the spectrum up to 100GHz and has 
its maximum at lower frequencies compared with the one- 
component models. The compact component II provides the 
sub-mm fluxes. This picture is well distinguished from the ra- 
dial non-uniform self-absorbed source discussed by de Bruyn 
(1976) which is optically thick at any wavelength, while our 
model is optically thin between 4 — 100 GHz. The drop of the 
electron energy inside the compact core is an essential property 
of the model and is accompanied by a significantly increased 
magnetic field in the core. The physical parameters are given 
in Tab. 0. 



7. Inverse-Compton Spectrum 

The first consequence of the interpretation of the radiospectrum 
as optically thin synchrotron radiation from relativistic elec- 
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Fig. 9. The luminosity of Sgr A* from the radio to the soft X- 
ray regime for an assumed distance of 8.5 kpc. The ROSAT 
flux at 1.85 keV is included. The two-component synchrotron 
spectrum from a hot (T = 6.0 10 11 K) extended envelope I 
and a cooler (T = 1.6 10 11 K) compact core II (the same as in 
Fig. ^ is fitted as a solid line. It includes SSC and absorption 
as discussed in Sec. [7]. The intrinsic spectrum of the core (en- 
velope) corresponds to the dashed (dotted) line. The physical 
parameters are given in Tab.0. 



trons is the synchrotron self-comptonization (SSC) process. In 
this process, only synchrotron photons are considered as the 
low energy radiation scattered by relativistic electrons. The 
electron energy distribution n e (pf) is also obtained from the 
radiospectrum. From arguments on total cross sections for in- 
verse Compton scattering (Blumenthal & Gould 1970), we can 
estimate the emitted power due to comptonized synchrotron 
photons to be 



Pi 



Comp 



£/ph p 
U B ync 



since synchrotron radiation can be considered as Compton scat- 
tering of the virtual photons of the static magnetic field. The 
emitted power is scaled by the relative energy densities of the 
photon field and the static magnetic field. Taking an intrinsic 
synchrotron luminosity of 465L Q for a gaussian electron dis- 
tribution and a source radius of 3.75 10 13 cm, we obtain an en- 
ergy density of £/ph = 3.4 10~ 3 erg cm~ 3 , which falls short of 
the energy density of the magnetic field by more than two or- 
ders of magnitude. Thus we get a Inverse-Compton luminosity 
of 



Po 



omp 



4.1 10 



"4.0G" 


2 


'3.7510 13 cm~ 


-2 

erg sec -1 


B 




R 





The spectrum of the scattered photons can be calculated us- 
ing the Thomson cross section <tt, since the synchrotron pho- 
tons are soft in the electron rest-frame jhu <C m e c 2 . For an 
isotropic electron and synchrotron-photon distribution (at least 
on average), the inverse-Compton flux is 



F v = ——o T chv 

l07T 



d7 n e (7) / dv s n s (v s )f(x) 



(22) 



Fig. 10. The luminosity of Sgr A* for assumed distance of 8.5 
kpc from the radio to the soft X-ray regime. The model spec- 
trum accounts for synchrontron and single inverse-compton 
emission. The ROSAT and Einstein flux is included. The 
dashed line is the spectrum including absorption (see text). 

where we have used the distribution f(x) (e.g Blumenthal & 
Gould 1970) for single scattering events, defined as 

f(x) = 2x ln(x) + x + 1 - 2x 2 x = — V —— 

and the index s refers to the synchrotron photon distribution 
and frequency. The mean density of incident photons in a ho- 
mogeneous spherical source of radius R with r = r(i?) = 
a v R and the source function S Vs is given by 



4tt 



1 1 



-2t 



n s (y s ) - .... ( , „ 
chv s \ t It* 

For most distributions, the integrals in Eq. ( |22[ ) can not be eval- 
uated analytically. The calculated spectrum shown in Fig. [l(] 
is produced by a gaussian electron distribution according to 
model IV of Tab.|l]. The measured flux density in the 0.8 — 2.5 
keV range reported by Predehl & Triimper (1994) for ROSAT- 
PSPC without correction for absorption is included in Fig. [Io| 
We take the upper limit for the ionized hydrogen column den- 
sity towards Sgr A* (B96) N H = 2.2 10 21 cm" 2 and obtain 
a cut-off due to photoelectric absorption (Morrison & McCam- 
mon 1983) at approximatly 1 keV. SSC and external absorption 
with the predicted column density from the turn-over at 1 GHz 
can account for the ROSAT measurement. Other measurments 
with Einstein (Watson et al. 1981) and ART-P (Pavlinsky et al. 
1994) are not consistent with this interpretation. The assumed 
TVh column density is much too low to account for the visual 
extinction of approximately 30 mag corresponing to a canonical 
hydrogen column density of = 6 10 21 cm~ 2 towards the 
galatic center. 



8. Summary 

Synchrotron radiation by quasi-monoenergetic electrons is 
considered as an explanation of the radio spectrum of the galac- 
tic center source Sgr A*. Therefore, we developed a detailed 
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treatment of the emission process of single electrons in a ho- 
mogeneous magnetic field to obtain a correct synchrotron spec- 
trum even for mildly relativistic electron energies. We showed 
that the emissivity of a single electron is modified due to the 
relativistic motion of the electrons with respect to the observer 
compared with the standard theory. This effect is removed 
when a stationary and isotropic distribution is considered as 
shown in Appendix B. The spectra of monoenergetic electrons 
are integrated for various electron distributions covering a thin 
shell in momentum space. These quasi-monoenergetic elec- 
trons in a homogeneous and isotropic source provide a fairly 
good fit to the time averaged radiospectrum of Sgr A*. The 
resulting self-compotonized X-rays are sufficient to explain 
the ROSAT-observations with a very low hydrogen column 
density. This hydrogen column density is consistent with the 
upper limit of ionized hydrogen in the vicinity of Sgr A* (B96) 
infered from free-free absorption. Additional material is re- 
quiered to account for the visual extinction. This may also lead 
to enhanced absorption in the soft X-ray range and the syn- 
chrotron self-comptonization would not be sufficient to explain 
the flux measured by ROSAT. 

Inverted radiospectra showing a synchrotron flux propor- 
tional to v - 3 are also reported for the centers of M81 (Reuter 
& Lesch 1996), NGC 1068 (Wittkowski et al. 1997) the archety- 
pal Seyfert 2 galaxy and M 104 (Jauch & Duschl, in prepa- 
ration). The spectra of all these sources can be interpreted 
as due to optically thin synchrotron radiation from quasi- 
monoenergetic electrons. This suggest a new common feature 
for the centers of normal and active galaxies. The differences 
in the radio regime arise from the attainable electron energies 
and the source size. 
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Appendix A: The Radiation Field at the Observer 

The electric field at the position of the observer can be ex- 
pressed in terms of source properties at the retarded time (Jack- 
son 1962) 



E = \s 



git) 



cR* 
n x (n 



/3)x/3 



(1 



n 



(Al) 



(A2) 



We take the viewpoint of an observer at an angle 9 between 
the line of sight and the magnetic field lines. Assuming a large 
distance between the source and the observer we approximate 
the distance at the retarded time t = t + R(t) /c as 
R{t) ps R- nr(r) . (A3) 
The vector 

^ c / sinacosi/) 
r(r) = — sin a cos 9 sine/) — tuj s sva9 cos a | (A4) 
8 \ sin a sin 9 sin d> + raj* cos 9 cos a 



describes the path of the spiraling electron with the azimuth an- 
gle (p = u> s t inside the source and n = (0, 0, 1) the line of sight 
in a coordinate system with the observer in the z-direction. If 
we take 9 to be a constant for sufficiently long times, it fol- 
lows that w(t) becomes a periodic function in r. This forces 
us to expand the electric field as was previously stated by Shu 
(1991). Thus the spectrum of a single electron decomposes into 
discrete emission lines with Fourier coefficients given by 



27 



dt g(t) exp[— imuit] 



(A5) 



and the time average of the observed power is provided by Par- 
seval's theorem 

r— 2 00 

P(fi) = v I " dt P(Sl, t) = ±- \9m\ 



o 



(A6) 



and the symmetry of g rn with respect to m — > — m. We trace 
back the integral to the retarded time and, following the deriva- 
tion given by Bekefi (1966), we obtain 

— imtot ~ —im{uj s T — ?/;sin(w s r)) (A7) 

for the argument of the exponential function. This provides us 
with a parameter reflecting the geometry, ip, and connects the 
frequency of the spiral with the frequency measured by the ob- 
server u). The derivative of the light cone conditon 

dR(r) > 



df = dr 1 



idr 



(A8) 



allows a substitution of the integration variable t by r in the 
same way as Jackson (1962 Cap. 14.5). Using approximation 
( A3 ) and the identity 

d n x [n x j3] 



dr l-nf> ^-^^ 
we perform a partial integration of 



) and get 



9m = 

with 



imoj 
2tt 



dr gi exp[— imuit] 



(A9) 



(A10) 



(All) 



sin a sm < 
<l = I sin n c:o> 9 c<>- < • • > 


t = (t(1 — f3 cos a cos9) — Puj~ x sin a sin#, sin(a; s T)) 

The contribution from the boundary in the partial integration 
vanishes if the argument is periodic with period 2iruj~ 1 which 
determines the frequency u> = (1 — /3cos# cosa)~ 1 ui s seen 
by the observer. Changing the variable of integration from r to 
4> we get 



2ire 



> <7i exp[— im{(j) — ip sin cj>)] 



(A12) 



with ip and £ defined in (|J). The argument of the exp-function 
in ( \12) is antisymmetric in <fr — > — </> and the ^-component of 



gi is also antisymmetric. The x-component of g m is given by 
the symmetric part of the integrand as 



P—t sma 
cos((m - 



d(f> [cos((m — 1)4> — mipshnj)) 



\)(j> — mip sin (f>)] 



(A13) 
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Taking the integral representation for Bessel functions of inte- 
ger order (e.g. Abramowitz & Stegun 1972 [9.1.21]) 

r 

7rJ m (z) = / &4> cos{z sin 4> — mcp) (A14) 
Jo 

and the recuiTence relation J m _\{z) — J m+ \{z) = 2J' m (z) the 
x-component of g m becomes 

sin a J' (rmp) (A15) 

The y-component of g\ is symmetric and the part contributing 
to the y-component of g m is 

.mw s 



2£ 



l P[2 cos a sin 9 J m (rmp) 



— smacos6(J m _i(mip) + J m+ i(mip))] (A16) 

With an additional recurrence relation J m -i(z) + J m +i(z) — 
2mz~ l J m (z) the y-component of g m can be rewritten as 

. moj s (3 cos a — cos 9 



sm t 



-J m {mip) 



(A17) 



The two representations ( A15 ) and (417) for the two compo- 
nents of the electric field are presented in Eq. (0). 



Appendix B: The Volume-Integral of Recieved Power 




Fig. Bl. Geometry of the the source region and choice of coor- 
dinates. 



B.l. Geometry and Coordinates 

We consider the emission of synchrotron radiation by rela- 
tivistically moving electrons in a homogeneous magnetic field 
which is inclined to the observer by an angle -d. The source 
is located far from the observer and is unresolved, so that 
the recieved power does not depend on the position of the 
emitting electron inside the source. Due to the motion of 
an individual electron along the magnetic field with velocity 
v\\ = cj3«, this electron is seen by the observer for a time 
At = ± | cos$) if x is the length of the elec- 

tron's path along the magnetic field during the observation. 
The time interval is shortened for electrons approaching the 



observer and prolonged for those which are receding. For cal- 
culating the observed power, Scheuer (1968) showed that the 
received power of an individual electron is reduced by a factor 
(1 ± /3ii cos$) since it contributes to the total power only for 
a velocity-dependent time interval. When dealing with a distri- 
bution function in phase space discribing the electrons in the 
source we show here that this argument must be based on the 
invariance of the distribution function to get the same result. In 
order to simplify the discussion we consider the source to be 
a box of length 2x$ with the magnetic field directed along the 
x— axis and the centre of the coordinate system centered in the 
box. Suppose we have two fluxes of electrons all of the same 
energy. One moving in the negative x— direction with velocity 
—dm and the other in the positive direction with v\\ . Synchrotron 
radiation should only be emitted while the electrons are inside 
the box and the distance between two neighbouring electrons 
moving in the same direction is x$ /N. The distribution of elec- 
trons can be written as 



f(x) = Q(x - \x\) 



n 

53*(s + m^-c0||t) 



(Bl) 



For each emitting electron it takes a time Ai e for the radiation 
to arrive at the position of the observer, given by 



c 1 (I — X cos 1?) 



(B2) 



which depends on the distance between source and observer I 
which is omitted in the following analysis. 



B.2. Normalizing the Distribution 

In general the distribution function and its integral in phase 
space is an invariant due to changes in the variables. This can 
be expressed as 



At 



At I A 3 r A 3 p f(r,p,t) = AN At 



with the knowledge of the energy equation E 2 = p 2 c 2 + m 2 c . 
We have AN electrons on average in the volume of interest and 
the complete distribution function of our problem reads 

/ = [5 3 (p-p )+S 3 ( P + Po )]e(x ~\x\) 



<% - Vo)S(z - zo) 



x — n- 



xo 



Px 
rn, 



(B3) 



where we used the momentum po = (vu m e , 0, 0). Performing 
the simple integrals gives us 



AN At = 



[ At 


1*1 


/At J 





Ap x f(x,p x ,t) 



If we want to express the integral in terms of the arrival time of 
signals at the observer t = t + - cos we have to change the 
integration variables, resulting in a factor 

d{x,p x ,i) _ px 
d(x,p x ,t) m e c 



■ cosz? 
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with x — x and p x — p x . So, whenever we calculate mean 
values of observables P(x,p x , t) in terms of x, t, we have to 
take 

P(t) = 



dx 



X 

c 



p_ 



(B4) 



(B5) 



N " v c 

m 

= 1 + P\\ cos-d P + = P(x, +p a ) 
= 1 - /3|| cost? P_ = P(af, -po) 

B.3. Mean observed synchrotron power 



If we know the received power for the two different types of 
electrons to be P + , P_ we can integrate to get the total received 
power. If we collect terms containing x in the argument of the 
(5-functions in (B4), we find 



P(f) 



dx 



P+Zi 



HZ+x 



x 
I — 
N 



■C0\\t 



(B6) 



for the total power recieved. To evaluate the spatial integral we 
transform the argument of the S— functions 



P = P4 



/XQ 




•?'o 



Z+N & 



x + m 



(B7) 



Z-N Z- J 

according to / S(f(x))dx = J2 n S ( x ~ x n)\\f {x n )\\~ l . The 
conditions under which the electrons, which are counted by n 
and to, contribute to the integrals is 

x J \x J 



xo 



Z- 



N < to < —8»t + Z- 
v x 



N 



An = 2£ + N & Am = 2£-N (B8) 
and we conclude that, for any instance of time, the number of 
electrons contributing to the integrals is 2N£ + and 2iV£_ re- 
spectively. In this way the argument given above for the differ- 
ent time intervals is recovered. We see that the total received 
power is changed by this counting argument only because the 
change in the time as a variable t — > t is compensated by the in- 
variance of the integral of the distribution function in T-space. 
This results in 

P = 2N(1 + /3|| cos tf)P+ +2N(1- /3|| cos i?)P_ (B9) 
The results can alternatively be derived by multipling the num- 
ber of observed electrons of each kind given in (B£) by the 
power recieved from these electrons. In this way, the result 
based on counting individual electrons (Scheuer 1968, Rybicky 
& Lightman 1979), can be recovered. 
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